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ABSTRACT 

The R Coronae Borealis (RCB) stars are hydrogen-deficient, variable stars 
that, along with the probably related H-deficient Carbon stars (HdCs) and Ex- 
treme Helium stars (EHes), are now believed to be the results of He-CO WD 
mergers. Their H-poor, C-rich atmospheres are characterized by very unusual 
isotopic ratios and elemental abundances, which provide stringent constraints on 
the physical processes during and after the double-degenerate merger. These in- 
clude extremely low oxygen isotopic ratios, 16 0/ 18 ~ 1 — 10, 12 C/ 13 C > 100, 
and enhancements up to 2.6 dex in F and in s-process elements from Zn to La, 
compared to solar. As shown before such O-isotopic ratios cannot result from the 
dynamic double-degenerate merger phase, and we therefore investigate now the 
role of the long-term post-merger evolution and nucleosynthesis. We construct 
post-merger ID spherical models based on realistic merger progenitor models and 
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the outcome of our previous hydro dynamic simulations and follow their evolution 
into the giant domain of the HR diagram where RCB stars are located. Along 
with convection zones that develop in the stellar envelope, we adopt a model for 
extra mixing that is meant to represent processes driven by rotation originating 
in the dynamical merger. We perform complete multi-zone post-processing nucle- 
osynthesis simulations for these stellar evolution models, using a nuclear network 
that includes over 1000 isotopes. Our models reproduce, for the first time, the 
full range of the observed abundances for almost all the elements measured in 
RCB stars: 16 0/ 18 ratios between 9 and 15, C-isotopic ratios above 100, and 
~ 1.4- 2.35 dex F enhancements, along with enrichments in s-process elements. 
The nucleosynthesis processes in our models constrain the length and tempera- 
ture in the dynamic merger shell-of-fire feature as well as the envelope mixing in 
the post-merger phase. Depending on model parameters s-process elements origi- 
nate either in the shell-of-fire merger feature or during the post-merger evolution, 
but the contribution from the AGB progenitors is negligible. The post-merger 
envelope mixing must eventually become inefficient at ~ 10 6 yr after the dynamic 
merging phase but before the star enters the RCB phase. 



Introduction 



RCB st ars are rare F-G type supergiants belonging to the old bulge population in 
the Galaxy fjCottrell fc Lawsonlll998l ; iTisserand et al.ll2008f ). Their atmospheres are mainly 
composed of He, almost entirely depleted in H yet rich i n C. They are n ear-solar mass single 
stars (except for DY Cen which is in a binary system, iRao et al.l |2012| ) that have single or 
multimode pulsation cycles. The pulsations are considered to play a role in the deep, irregular 
brightness declines of up to 8 mag nitudes jCkytoiJll996l ). These declines are attributed to 
clouds of dust ejected fro m their atmosp heres, the evidence for which is an omnipresent IR 
signature around the star (jClaytonlll996l ). Such deep declines and IR excesses have not been 
observed in the otherwise chemically similar Hydrogen-deficient Carbon (HdC) stars or the 
hotter Extreme Helium (EHe) stars. 

Other than the declines, for which they were first noticed, RCB stars have been an 
enigma because of their anomalous chemical signatures and the possible evolutionary path 
(or paths) that may have led to their current state. The surface layers of RCBs are pre- 
dominantly made of helium (mass fraction Xh r =0.9 8) and are carbon rich with C>N>0. 
They also have 12 C/ 13 C ratios larger than 40-100 (IHema et al.l 120121 ) (compare with the 
CN-cycle equilibrium val ue ~ 3.5 that is expe cted for H-burning ashes), F enhancements in 
the range of 1.0 -2.7 dex (jPandey et al.ll2008l ) and considerable overabundances of s-process 
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elements compared to solar values (jAsplund et al.l 120001 ). This peculiar surface composition 
is indicative of material processed by H-burning and (partial) He-burning. The 16 0/ 18 
ratios in HdC and RCB stars have been a focus of many recent studies, sta rting with the 



disco very of an extremely low value of 0.83±0.02 in the HdC star HD 137613 (jClayton et al. 



20051 ). as compare d to the solar value of 500. RCB stars, as a class , have 16 0/ 18 ratios 
between 1 and 20 (jClayton et al.l 120071 ; iGarcia-Hernandez et al.l 120101 ) . probably the lowest 
ones measured anywhere in the Universe. The goal of this work is to construct post-merger 
models that can reproduce the unique chemical signatures of RCB stars, particularly their 
low 16 0/ 18 ratios and F enhancements. 

The EHe stars have similar compositions (barring a few exceptions) as RCB stars, 
suggesting a common evolutionary origin and link. The surfaces of EHe stars are too hot to 
show CO bands, therefore their 18 0/ 16 ratios cannot be measured. As one of our primary 
goals in this paper is to model the low 16 0/ 18 ratios seen in RCB stars, we exclude EHe stars 
from our discussion. It should be noted however tha t EHes do have high F enhancements 
reported, like those in RCB stars (jPandey et al.ll2008f ). 



65 RCB s tars have been found in the Milky Way and 25 in the Magellanic Clouds 
(jClaytonll2012l ). Two competing scenarios have been put forth to explain how RCB stars 



may have formed - the Final He- Flash (FF) model dRenzini 



(DD) white dwarf merger model (IWebbinklll984j ; llben et al. 



1990h and the Double Degenerate 



1996J). In FF model, a post-AGB 



star experiences a He-shell flash during the pre-white dwarf phase and evolves back into a 
giant configuration. In the DD merger scenario, a CO white dwarf and a less massive He 
white dwarf (WD) in a close binary system, formed as the result of one or more common 
envelope episodes, finally merge over a dynamic period of ~ 100 - 1000 s. In terms of number 
statistics, the expected number of RCB stars resulting from the merger of CO and He 
WDs is very clos e to the estim ated 5700 RCB stars obtained by extrapolating the LMC 
RCB population (jClayton! I2012T ). which implies that only a small fraction has been found 
so far. Computational models and observational evidence also lean toward the DD model 



being the more plausible path out of the two scenarios to for m RCB stars (jClayton et al. 



20071 ; lLoren-Aguilar et al.l 120091 ; I Garcia- Hernandez et al.l |2009|) . The main nucleosynthesis 
evidence comes from the 12 C/ 13 C and 16 0/ 18 ratios in RCB stars. 

According to a single- zone analysis, rapid mass-transfer events in He+C O WD binaries 



could lead to thermodynamic conditions in which such ratios would originate (jClayton et al. 



20071 ). Specifically, 18 can be produced by partial He-burning via 14 N(a, 7) 18 F(/3 + ) lo 0, 
without being destroyed by a further a-capture to 22 Ne for a temperature in a range of 
1.2-1.9xl0 8 K at a density of 10 3 g/cm 3 and on time scales of ^250-320yr. At these 



,18, 



temperatures, 13 C is consumed by He-burning to form 16 via 13 C(a, n) 16 0, which leads 
to the high C-isotopic ratios observed in RCB stars, while the neutrons released from this 
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reaction could produce s-process elements. 



The thermal properties of a final-flash event are equivalent to those of an AGB ther- 
mal pulse, with temperatures in the cony ective He-burning shell as high as 3 x 10 s K 
(IWerner fc Herwigll2006l ; iHerwig et al.ll201ll and ref. therein) maintained for a long enough 



time to process 18 and 19 F to 22 Ne (t he latter via the reaction 19 F(a, p) Ne ). This may ex- 



plain the lack of 18 O (from CO bands, fevres et alill998l : beballe et alil2002h ) and F I bands 
( Pandey et al. 20081 ) in the spectrum of Sakurai's object (V4334 S gr) which is believed to 
be a final-flash star of the very late thermal pulse (VLTP) variant (jHerwigll200ll . and refer- 
ences therein). In the VLTP event, 13 C is produced by proton inge stion into the convection 
zone and their capt ure by 12 C resulting in low 12 C/ 13 C ratios of ~ 5 (IHerwig fc Langer 112001 



(lAsplund 



1997 



Herwig et al.ll201lh. w hich compares well with the observed value of 1.5— 5 in Sakurai's object 



1999). VLTP star models also predict Li enrichments ( IHerwig &: Langer 



20011 ). as observed in Sakurai's object, but there are only four Li-rich RCB stars. 



A compact system of two WDs evolves from a pair of main sequence stars that are close 
enough to interact with each other through the exchange of gravitational wave radiation and 
magnetic stellar winds, causing a loss of orbital angular momentum from the system and thus 
reducing the orbital separation be tween the stars. Such c lose d ouble-degenerate WD systems 
form over a period of ~ 10 9 years (jlben fc Tutukov!ll984l . 119871 ) through one or more episodes 
of mass-transfer (including at least one common-envelope interaction). The final merging of 
the two WDs occurs on a dynamic timescale (100-1000 seconds), with the less massive WD 
(the He WD) being tidally disrupted by the more massive one (the CO WD). Smooth Particle 



Hydrodynamic (SPH) simulations t hat follow this merging phase (lLoren-Aguilar et al.ll2009 



Raskin et al.ll2012l ; IStaff et al.l 120121 ) show that the He WD material rapidly inf ailing on the 
CO WD causes the surface of the latter to be shock-heated, resulting in a hot corona feature, 
while the rest of the He- WD mass settles in a disk-like structure that extends to large radii. 

One-dimensional He-rich homogeneous stellar evolution models with C, O and trace 
amounts of H showed promising evolutionar y tracks into the RCB r egion in the HR diagram 
(IWeisd 119871 ). In an alternative approach, ISaio fc Jeffery I (120021 ) represented the merger 
phase through accreting He- WD material (of a pure helium composition with a small H- 
rich envelope mass) onto a CO WD at the rate of lO -5 M /year. Although the location of 
the final model in the HR diagram was consistent with that of EHes and RCB stars, the 
accretion rate wa s much smaller than the exp ected rapid mass-transfer rate of ~ lO 3 M /year 
(jStaff et al.ll2012l ; lLoren-Aguilar et al.l 120091 ) during the final merging phase of WDs. 



The SPH simulations of the dynamic merger phase of He and CO WDs , which were 



performed over a range of mass combinations by lLoren-Aguilar et al.l ( 120091 ). included a 
14-species nuclear network. The abundance yields in the hot coro na of the final models 



of masses > 1.2 M showed enrichments of Ca, Mg, S, Si, and Fe. lLongland et al.l (120111 ) 
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carried out a one-zone post-processing nucleosynthesis analysis using the temperature and 
density conditions from the corona of the final model. The AGB model used to build their 
initial nuclear composition contained 7 nuclear species from 1 H to 22 Ne. They reported that, 
depending on the depth of mixing within the corona, the 16 0/ 18 ratio had decreased to 
~ 19 -370. However, their total C and O abundances exceeded the upper limit for RCB 
stars and their F predictions were much lower than the observed values, while predictions 
for s-process element production were not available. 



Jeffery et al.l (120111 ) have investigated RCB surface abundances by using a cold mixing 



recipe, i.e. without nucleosynthesis during either the hot merger phase or the post-merger 
evolution. A 1-D 'cold merger' model involved mixing different proportions of masses of He 
and CO WDs, and a 'hot merger' model was built by inc orporating the resul t s obta ined 
in the SPH simulations of He-CO WD merger systems by lLoren-Aguilar et al.l (120091 ) . In 
both cases, the elemental abundances obtained from their mixing recipe, barring those of 
C, N and O (and S in the hot merger model), fell short of the observed ones in RCB stars, 
particularly the s-process elements, which were assumed to originate exclusively from the 
progenitor AGB star. They have concluded that even additional nucleosynthesis that may 
occur during the actual merger of He and CO WDs cannot help to reproduce the observed 
abundances of RCB stars. 



Continuing the work by IClayton et al.l (120051 ) and IClayton et al.l ( 120071 ). IStaff et al. 



(120 12[ ) (paper I) presented grid-based hydrodynamic merger simulations for a range of mass 
ratios (q) of CO a nd He WDs (of the total m ass 0.9 M ). A Shell of Fire (SOF) (similar to 
the hot corona of lLoren-Aguilar et al.l 120091 ) was identified in the merged object, and that 
feature was analyzed for the potential formation of low 16 0/ 18 ratio material. Single- zone 
nucleosynthesis calculations using the T, p conditions of the SOF of three cases were made 
for a range of initial abundance distributions derived from a variety of appropriate AGB 
progenitor models. Those were taken from the NuGric0 library of detailed post-processed, 
fully evolved models of solar metallicity that ranged from the early AGB to the advanced 
thermal pulse stage (Set 1.2, Pignatari et al., in prep). The composition of the He WD was 
taken to be uniformly mixed and w as extracted from a low mass RGB star. In contrast to 
the results of lLongland et al.l (120111 ). we found that within the dynamic phase, the 16 0/ 18 
ratio in the SOF was of the order of 1000. The difference in the 16 0/ 18 ratios between 
the two simulations may arise from the difference in the total mass considered. In any case, 
according to our results, the low observed 16 0/ 18 ratios of RCB stars cannot originate from 
the dynamic merger phase, but if the temperature and density conditions were maintained 
for a longer period of time, then 16 0/ 18 ratios between 4 (q = 0.7, 10 9 s) to 30 (q = 0.5, 
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10 6 s) could be obtained. 



Following the results of paper I, we now explore the the long term post-merger evolution 
with the goal of identifying conditions under which this evolutionary phase could reconcile 
the present stark discrepancy between models and observed surface abundance of RCB stars. 
We are motivated by the fact that merging WDs are considered to be paths to the creation 
of many other interesting stellar objects, such as Type l a supernov ae and neutron stars 
jHanlll998 : lNelemans et al. IhoOlatlNapiwotzki et al. IhoOlh . sdB stars jSaio fcJeffervlhoOol: 
Han et alj[2002[ ). AM CVns JSolheiml 1201^ ) and even planetary nebulae (IDe Marco et al.l 



20091 ). The rich and detailed observational picture available for RCB and HdC stars gives an 
opportunity to constrain DD merger and post-merger physics in general, thereby providing 
a quantitative understanding of the objects mentioned above. 

The long term evolution of such merged sy stems may be characterized by the formatio n 
of an accretion disk from the secondary WD (jYoon et al.l 120071 ; Ivan Kerkwijk et al.l |2010| ) . 
Alternatively, the post-merger object may undergo a more yiscous, star-like evolution with 
the K eplerian rotation in the envelope eventually dying out (jShen et al.ll2012c ISchwab et al. 
20121 ). We follow the latter scenario for the post-merger evolution, and we assume the merged 
object to evolve as a hybrid C/O-He star undergoing rotation-driven mixing in its envelope. 
We construct spherically symmetric one-dimensional models of merged WDs, with a range 
of initial conditions and follow their evolution into the domain of the HR diagram where 
RCB stars are observed to lie. Along with convective mixing, these models also include 
an induced, continuous mixing profile assumed to result from the rotation of the merger- 
remnant. We then perform a multi-zone post-processing nucleosynthesis analysis of these 
models and compare their surface abundances with those of RCB stars. 

This paper is organized as follows. Section 2 describes our simulation methods, where the 
progenitor evolution of the system, the construction of the initial models, a brief description 
of the codes we used, the formulation of the mixing profile, and the four post-merger evolution 
cases that we studied are presented. This is followed by Section 3, where the stellar and 
nucleosynthesis processes of the model during its evolution are described and compared with 
each other and with the RCB star observations. We close with the discussion and conclusions 
in Section 4. 



2. Methods 

The post-merger simulations are building on our previous dynamic merger simulations 
(paper I). They involve the construction of appropriate initial models for the post-merger 
evolution, starting with the abundance profiles that take into account information for all 
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observed species from realistic progenitor models as well as our best estimates of nucleosyn- 
thesis in the dynamic merger phase (i.e. in the SOF). The initial abundance profiles are 
homologously imposed on a homogeneous zero-age main-sequence models of the same mass, 
so that within a thermal time scale the stellar structure is determined by and consistent with 
the post-merger abundance profile. The stellar evolution of these objects is followed into the 
RCB region of the HRD, and then post-processed with the same complete nuclear network 
(using the NuGrid PPN codes) that has been used to generate the progenitor models. In this 
way our models take into account nuclear production during the progenitor evolution, the 
hot merger phase (paper I) as well as any production during the post-merger phase, which 
- as we show here - is essential to reproduce the observed abundances of RCB stars. 



2.1. Progenitor evolution and outcome of dynamic merger phase 



The progenitor evolution of a double-degenerate merger evolves from a close binary 
of low- and intermediate-mass main sequence sta rs, and involves a common envelope (CE) 
event when the primary becomes a giant star (e.g. Ilben et al.lll996l ). CE eve nts are expected 



to have occurred in clos e binary systems with at least one WD component ( iDe Marco et al. 



20 111 ; iPassy et al.ll2012l ). The outcome of the second CE event is a close binary with CO 
and He WDs. This progenitor evolution matters here only in as much as it determines 
which RGB and AGB models should be considered when constructin g the initial abun dance 
distribution out of their combined He and CO cores, as described in [Staff et al.l (120 121 ). 



The final merging of the two WDs takes place within a dynamic timescale of 10 2 -10 3 s 
and results in a structure consistin g of a CO core surrounded by a shock-heated layer and 



an en velope in Keplerian rotation ( jLoren-Aguilar et al.ll2009l ; iDan et al.l 120121 ; iRaskin et al. 



20121 ). In our hydrodynamic simulations reported in paper I, this hot and very dense SOF 
lasted until the end of the simulations. The SOF is the region where nucleosynthesis can 
occur within the dynamic time scale. The final nuclear production of the SOF will depend 
on how long it lasts beyond the end of simulations. 

The post-merger evolution is determined by the interplay of thermal dissipation, angular 
momentum transport induced by rotational instabilities and magnetic fields, and the nuclear 
time scale of H- and He-burning shells. How these ingredients interact is still uncertain. One 
scenario involves accretion of the post-merger d isk remnant onto the core over a timescale of 
either 10 5 years (jYoon et al.ll2007| ) or few hours ( Ivan Kerkwijk et al.ll2010l ). Alternatively, the 
post-merger object couples the rigid primary WD and the tidally disrupted secondary WD 
undergoing Keplerian rotation through viscous mechanisms (from magnetohydrodynamic 
instabilities) that stimulate t ransport of angular momentum over a period between 10 4 and 
10 8 seconds (jShen et al.ll2012l ). These shearing forces may also increase the peak temperature 
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of the SOF/cor ona of the remnant b y a factor of ~2 and work on expanding the material 



to larger radii. ISchwab et al.l (120121 ) show that, towards the end of such a viscous phase, 
the remnant evolves towards a spherically symmetric, almost shear-free steady state with a 
thermally-supported envelope in solid body rotation. These authors expect the subsequent 
evolution to be star-like, growing into a giant star, with the luminosity from its internal 
nuclear-burning driving convection in its cool extended envelope. Such evolution represents 
a path from the merger of white dwarfs to the formation of RCB giant stars. 

Following the latter scenario, we construct the initial composition profile for our long- 
term post-merger simulations to be comprised of a CO core, an SOF, and an envelope with 
rotationally induced mixing, initially representing the former He WD. 



2.1.1. The post-merger composition profile 

The first step in building the internal composition of the post-merged object is to choose 
the progenitor AGB and RGB stars of the CO and He WDs respectively. These progenitor 
models determine when the CE interactions will take place and therefore cannot be uniquely 
determined a posteriory. We consider two out of the three CO- WD progenitor models intro- 
duced in paper I. Accordingly, the CO WD is the progeny of either early-AGB (E-AGB) or 
advanced-thermal pulse AGB (A-TP AGB) star. The difference between these progenitors 
is in the amount of s-process elements in the outer layers of the CO WD@, as well as the 
thickness of the He-rich layer. In both cases, the CO WD has a H-rich envelope of the 
order of 1O _4 M . The RGB model for the He WD has a more massive H-rich envelope of 
~ 10~ 3 M Q . The masses of the different layers of the WDs are listed in Table [TJ 



2 However we find that this s process material in the progenitor AGB is not the origin of the observed 
s-process overabundances in RCB stars, cf. SectiorO 

Table 1. The WD masses used for the initial models of the post-merger evolution. 
M WD _ core is the mass of the WD core (Fig. [[]), while M env is the mass of the H/He- or 
H-envelope for the CO and He WDs, respectively. 



WD 


Progenitor (M init /M Q ) 


M W d/ M q 


AfwD-core/ Mq 


Af onv /lO- 4 M 


CO 


E-AGB (3.00) 


0.58148 


0.4876 (He-free) 


2.4 


CO 


A-TP AGB (2.00) 


0.61243 


0.5789 (He-free) 


1.7 


He 


RGB (1.65) 


0.3024 


0.296 (H-free) 


64.3 
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The various components of the pre-merger WD composition are combined by taking into 
account the dynamic mixing and the SOF nuclear burning, both derived from the merger 
simulations of paper I. This leads to a four-zone model consisting of the relatively cold, 
degenerate core of the CO WD, a thin buffer zone, the SOF itself, and a relatively cold 
envelope (Fig. [p. 

For each mass ratio, q = M He /M C o, a certain fraction of mass (/drc) is dredged up from 
the CO WD (Mco) and mixed with the He WD (Mn e ) material in the envelope. The amount 
of dredge- up defines the outer boundary of the core, m core (Eqs. [p. As this dredged- up mass 
fro m the He-rich lay ers of the CO WD mass (m dre ) is pulled up to the envelope (see Fig. 11 



of IStaff et al.l |2012| ). an equal mass of He WD makes its way into the CO WD. Thus a 
partially mixed zone is formed between the two WDs. A major portion of the mass in this 
zone is occupied by the SOF (msoF), while the rest is in the buffer zone between the core 
and the SOF (mbuffer)- A small fraction of the mass from the He WD (that penetrates 
the CO WD and forms the partial mixing_zone) resides in the SOF (mh e - S of), while a larger 
part finds itself located in the buffer zon 
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The initial SOF composition is thus dominated by the mass dredged up from the CO 
WD (m co _ so f). As in paper I, we assume that the exchange of masses during the dynamic 
merger phase occurs before the onset of burning in any given layer. Correspondingly, the 
SOF layer is subjected to nuclear processing according to the conditions in the SOF derived 
from the merger calculations (see Section f2. 1.21 for two cases considered for SOF processing; 
the temperature in the buffer layer of T fh 5 x 10 7 K is too low for any nuclear reactions 
during the dynamic phase). The result of this nuclear processing is then entered into the 
post-merger abundance profile. Finally, the envelope (m env ) contains the fraction of m dTe 
from the CO WD that is not present in the SOF and the remaining He WD mass. 



m dle = /drc x Mco 

m core = M co - m dre 

^buffer = m dre ~ "^he-sof (1) 
ITISOF = ^he-sof + "Ico-sof 

m env = M He - m co _ sof 



3 The detailed properties of this buffer zone are difficult to extract from the 3D merger simulations for 
the purpose of mapping them into the ID profile. However, it turns out that these details are not important 
for the final RCB abundance predictions. 
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Fig. 1. — Four-zone post-merger abundance profiles of the q =0.7, E-AGB CO WD with a 
dredge- up depth of 15% (case 1). 
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Motivated by the merger simulations, we assume that the He-rich envelope is completely 
mixed during the merger process. The most important consequence of this assumption is 
the uniform distribution in the post-merger envelope of H from the former H-rich He- WD 
envelope. After extracting the abundances from each zone of the progenitor models, the 
isotopic abundances of that zone are averaged over its mass. 



2.1.2. The four post-merger cases 

The key parameters of the post-merger models are the pre-merger WD mass ratio, the 
amount of dynamic dredge-up from the CO WD, the type of CO- WD progenitor model, and 
the temperature and duration of the SOF nuclear burning in the dynamic phase. We focus 
here on four representative combinations of these parameters. 

The dynamic merger simulations predict dredge-up depths ranging from 15% to 18% of 
mass into the CO WD for the three low-g cases. In the outer layers of the He-free core of the 
CO WD 16 O is the most abundant isotope. If material is dredged up from these layers the 
star becomes more enhanced in oxygen than in carbon at its surface, which is not observed 
in RCB stars. The other problem of dredging up high amounts of oxygen is that the star 
would have much higher 16 0/ 18 ratios than observed in RCB stars, even when post-merger 
nucleosynthesis effects are taken into account (see below). Therefore, in order to avoid an 
excessive surface enrichment in 16 O, the adopted mass of dredge- up from the chosen CO 
WD models is lower than the values obtained in the hydrodynamic simulations. We thus 
restrict the dredge-up depth in our models to layers just above the He-free core of the CO 
WD (ranging from ~5.4% in the A-TP to ~ 15% in the E-AGB CO WD model, Table EJ. 
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Table 2. Summary of the physical parameters of the four cases. These are the mass of the CO WD (Mco)> He WD 
(MHe), core (mcore), buffer zone (mbuffer), SOF (msoF) and envelope (m env ). Columns 4 and 5 specify the progenitor 
AGB models used to construct the CO WD and the percentage of dredge-up from within the CO WD. Column 10 
(tsof) gives the duration of burning of the SOF, while 11 and 12 indicate the period in the RCB phase and the colour 

representation of the cases in Fig. [lOj 



CASE 


Mco/ M Q 


M He /M 


AGB 


fdre (%) 


mCorc/ Mq 


TObuffcr/ Mq 


tosof/ M Q 


m cnv / M Q 


TSOF 


RCB age (10 5 years) 


colour 


1 


0.53 


0.37 


E-AGB 


15.0 


0.4505 


0.0295 


0.1000 


0.32 


short 


1.02 


orange 


2 


0.53 


0.37 


E-AGB 


15.0 


0.4505 


0.0295 


0.1000 


0.32 


long 


1.11 


red 


3 


0.53 


0.37 


A-TP AGB 


5.4 


0.5014 


0.0107 


0.0346 


0.352 


short 


0.97 


blue 


4 


0.60 


0.30 


E-AGB 


15.0 


0.5100 


0.0141 


0.1060 


0.2626 


long 


2.75 


green 



We consider two mass ratios from paper I, i.e. the q = 0.7 case (M o = 0.53 M and 
M He = 0.37 M Q ) with an SOF of T = 1.23 x 10 8 K and p = 5.02 x 10 4 gem" 3 , and the 
q = 0.5 case (M co = 0.60 M and M He = 0.30 M ) with an SOF of T = 2.42 x 10 8 K and 
p = 3.16 x 10 4 gcm~ 3 . The buffer zone has a temperature of approximately 5 x 10 7 K and 
the same density as the SOF, which is not hot enough for nucleosynthesis on the dynamic 
time scale. Similarly, the core and envelope are too cold for nuclear reactions during the 
merger. The duration of the high-T conditions in the SOF may extend beyond the time 
period covered by the hydrodynamic simulations, at the end of which the SOF is still hot. In 
addition, as pointed out in paper I, those simulations ignore any feedback of nuclear energy 
production during the hydrodynamic phase, which may affect the conditions and duration 
of the SOF. In order to project the impact of different SOF assumptions, we consider a short 
SOF with a duration of 8.2 x 10 5 s and a long SOF with a durat i on of 4.7 x 10 6 s, both of 



which are within the viscous timescale determined by IShen et al.l (120121 ) 



The four initial abundance profiles combine two CO- WD progenitor models, two mass 
ratios (q = 0.5 and 0.7) with their respective SOF parameters (T, p and the proportion of CO 
and He WDs in them), as well as the short or long duration of SOF burning, as summarized 
in Tabled 



2.2. Stellar evolution and nucleosynthesis simulations 

For the post-merger e volution, we use the version 3851 of the stellar evolution code 



MESA (IPaxton et al.ll201ll ). Initial models are created from a homogeneous 0.90 M (the 
total mass of the WDs in the hydrodynamic simulations of paper I) zero-age main-sequence 
model. This model is then relaxed to each of the four-zone post-merger abundance profiles. 

The star is evolved with a mass loss rate obtained by setting the coefficient of Blocker's 
wind formula to r] = 0.05, similar to what is used in the NuGrid AGB simulations. From 
the options available in MESA, we use the OPAL type I opacities for a solar abundance 
distribution. We have performed test runs with the more appropriate type II opacities 
that allow C and O enhancements from He-burning in the composition. In the end, such 
models are quite similar to type I opacity runs. However, models with type II opacities 
turn out to be much more difficult to converge, probably reflecting some of the instabilities 
that do lead to the irregular variability seen in RCB stars. Those test runs with type 
II opacities evolved to lower temperatures than their type I counterparts, but did not show 
differences in the nucleosynthesis aspects. In particular, the type II opacity runs did not show 
differences in the envelope convection properties that would have altered our assumptions for 
the adopted mixing profile (section I2.3p . Thus, the surface abundance predictions remain 
the same in models with and without CO enhanced opacities, and we choose to present here 
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a homogeneous grid of cases calculated with type I opacities. 



Once the stellar evolution calculations are completed, the track is post-processed for 
a complet e nucleosynthesis a nalysis. We use the NuGrid multi-zone nucleosynthesis code 
MPPNP ( IHerwig et al.ll2008l ) that processes each of the « 1300 zones for each of the w 3000 
time steps of the stellar evolution track with an adaptive nuclear network that includes all 
stable as well as all relevant unstable species (over 1000 isotopes) and their corresponding 
nuclear reactions. The post-processing code also mixes species after each time step, according 
to the mixing processes considered in the stellar evolution simulations. 



2.3. The mixing model 

RCB stars show abundance signatures at their surfaces that originate either in the 
SOF or in the nuclear processes during the post-merger evolution. Mixing in the envelope 
reaching into the nuclear processing region needs to bring such nuclear processed material 
to the surface. Early on in our investigation it became clear that convection alone would not 
be sufficient to provide this essential mixing ingredient in the envelope. 

The physical mechanisms of non-convective mixing processes that may occur during 
or after a merger are quantitatively not understood. In any case, the merger remnant is 



rotati ng, and it is reasonable to associate this rotation with mixing processes. IShen et al. 



(120121 ) estimate that differential rotation in such an object would last for at least 10 — 



10 8 seconds after the merging phase, after which it would evolve towards solid-body rotation 



(ISchwab et al.ll2012l ). The exact timescale for evolution towards solid-body rotation may be 
uncertain, but some mixing processes, such as meridional circulation, can be driven by solid 
body rotation as well. 

Considering the significant uncertainties regarding both the realistic rotation profile and 
the physics of rot at ion- driven mixing, we adopt an empirical approach to additional mixing 
processes in this study. This means that our primary goal is to identify a reasonable mixing 
law (a diffusion coefficient profile) that is capable of reproducing the observed abundance 
patterns in RCB stars. 

The diffusion coefficient in our mixing model (-D to tai) consists of convective mixing 
(Dmlt) predicted by the mixing length theory (MLT) and an additional mixing (-D a dd) 
component. The mixing is always restricted to the region above m m ; n ( = mcore in each case) 
(Eqs. [T]and Table [2]). -D a dd is the Eulerian representation of a constant mass diffusion co- 
efficient D Lag that drops sharply at the point where an entropy barrier appears due to the 
peak of nuclear energy generation. The location of this energy peak can be approximated 
as the mass-coordinate corresponding to a specified fraction of the 14 N abundance that has 



14 




0.40 
0.35 
0.30 
0.25 
t> 0.20 
0.15 
0.10 
0.05 



I I I I I I •! I I I I I I I I I I •] | | | 



logS 



0.00 



0.4 



0.5 



'•mi • i i i i i i i i i i 



--/ 



t 

'••••••< ,, < 



9.4 
9.2 
9.0 



0.6 0.7 

M r /M© 



0.8 



0.9 



- 8.8 


(-0 




OC 


- 8.6 




- 8.4 




- 8.2 




- 8.0 





Fig. 2. — The diffusion profile -D total at the timestep (iv) (Fig. |6]) for case 1 (black dot in 
Fig. E]), along with the MLT diffusion coefficient -Dmlt, in comparison with relevant physical 
parameters, such as the 14 N abundance (Xi4 Nk , panel 1), temperature (log T) and nuclear 
energy generation rate (e nuc ) (panel 2), as well as the radiative (V r ) and adiabatic (V ) 
temperature gradients and entropy (log S) (panel 3). 
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Fig. 3. — The diffusion coefficient profile (log -D total ) at the timestep (iv) (Fig. [6]) for case 1 
(black dot in Fig. [5]) along with the thermal diffusivity (log i/p). 
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still remained at this depth during He-burning, k = X( 14 Nk)/X( 14 N max ), where X( 14 Nk) is 
the mass fraction of 14 N at the mass coordinate where the mixing drops, and X( 14 N max ) 
is the the maximum 14 N mass fraction in the model at that time (Eqs. [2J Fig. [5]). 

The sharp drop of Aotai at the mass-coordinate is modeled by an exponential de- 
crease of the mixing efficiency. This transition from the extra-mixing zone to no-mixing is 
formally similar to how the depth-dependent convective overshooting is implemented in the 
MESA code, and like in the latter case, the width of this transition region is determined by 
the parameter / in the exponent of Eq. [2c] (see the steep ramp in Aotai i n Fig. [2). 



Aotai = m < m min (2a) 

Aotai = Add + £>mlt m >= m min (2b) 

n Aag /— (T — T k )) 

Add = j. — 9v? x ex P( 7^ ) m mhi <m<m k (2c) 

Add = Ta m > m k , (2d) 

where Tk is the temperature at the mass coordinate m^. This mixing model is specified 
with the following parameters: Aag = 4.5 x 10 51 g 2 /s, fc=0.003 and /=0.05. The resulting 
Eulerian diffusion profile is shown in Fig.[2j In Fig.EJ Dtotai is compared with the thermal 
diffusivity z/ T for one of our considered cases. It is seen that D total <C z/ T everywhere inside 
the star, which means that the additional mixing could be driven by a secular instability 
rather than by a dynamical one. It is interesting that the ratio between D to tai and u^, 



Dtotai ~ 10 2 z/t, is close to that used by iDenissenkov &: Pinsonneaultl (120081 ) to model extra 
mixing in low-mass red giants. 

The effects of varying the free parameters in Eqs. [2] are briefly summarized here. When 
we increase DL ag by more than two orders of magnitude, the mixing becomes so fast that it 
causes an excessive increase in 12 C and 16 abundances at the surface and a proportional 
decrease in 14 N that disagree with observations. There is also a risk for the star to become 
more enhanced in oxygen than in carbon, the opposite of which is observed in almost all 
RCBs. Making DL ag lower by a factor of 10 causes the mixing to become far too slow in 
dredging up 12 C and 16 to the surface. As the star enters the RCB phase, it becomes more 
N-rich than C, which again is not observed for any RCB star (except for one). 

The surface chemical abundances of the star vary within an order of magnitude when the 
value of the parameter k is varied between 10~ 3 and 5 x 10~ 2 . Below the mass coordinate 
m k , the diffusion profile drops very steeply. Allowing mixing to proceed deeper into and 
through the He-burning shell would cause an excessive dredge-up of 16 0. In addition, 19 F 
and 18 produced by partial He-burning would be destroyed in this case. On the other hand, 
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limiting the mixing to layers entirely above the He-burning shell would prevent these two 
species from being formed. Note that the location of diminishing mixing efficiency is 
physically motivated to coincide with the location of a steepening positive entropy gradient 
associated with the He-burning shell. 

A small ramp appears in -D tota i just as it nears m^. The beginning of the slope of the 
ramp and its steepness are decided by the rather sensitive parameter / in the exponent 
function of T. The ramp feature affects the content of 19 F, 18 and 16 at the surface (more 
discussed in section I3.2.2[) . The optimal abundances in the star at a given value of D^ ag and 
k are maintained when / is chosen between 0.03 and 0.06. 

Finally, we see indications from comparison with observed RCB star abundance features 
that the mixing law should be time dependent. In order to reproduce all the features, mixing 
must be strongly reduced shortly before the star enters the RCB region of the HR diagram, 
as will be discussed in section 13.3.11 

3. Results 

In this section we describe our results, starting with the nucleosynthesis in the shell of 
fire in order to define the initial conditions for the post-merger stellar evolution simulations, 
and finally we describe the nucleosynthesis in these post-merger evolution models. 

3.1. The Shell of Fire (SOF) 

The formation and evolution of the shell of fire feature in the dynamic merger simula- 
tions has been described in detail in paper I, where we have also presented nucleosynthesis 
simulations. Here we just need to briefly review and update those results with a focus on 
the conditions for the production of n-capture elements, such as the s-process elements. 

The temperature in the SOF depends on the mass ratio q of the merging white dwarfs. 
For the post-merger simulations presented here we consider two SOF cases: case 4 is based 
on a high-T (2.42 x 10 8 K, q = 0.5) SOF while the other cases are based on a low-T SOF 
(1.23 x 10 8 K, q = 0.7). 

In the low-T SOF 13 C forms early on at a level of ~ 10% by mass and exceeds the 
abundance of 12 C while 14 N reaches ~ 1%. In this situation neutrons are steadily released via 
the 13 C(a, n) 16 reaction with a neutron density reaching a peak value of ~ 3.16 x 10 11 cm -3 
(Fig.HJ). An interesting neutron- recycling effect is at play here, in which protons released 
by 14 N(n, p) 14 C are mostly creating another 13 C and then again a neutron. This explains 
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Fig. 4. — The evolution of the abundances by number fraction of 13 C(/100), Y, Ba and Ti 
during the burning of the SOF, along with that of the neutron density, in the low-T and the 
high-T SOFs. The vertical dashed lines indicate the short-SOF (8.2 x 10 5 s) and long-SOF 
(4.7 x 10 6 s) durations. The abundance of 13 C of case 4 (blue) is beneath the lower limit of 
this plot. 

why the 13 C abundance shown in Fig. @] decreases slower than expected considering just the 
a-capture reaction at that temperature. 

In any case, sufficient neutrons are available to cause a steady increase in the abundance 
of s-process elements along with nuclei such as Ti. This production of n-capture elements 
becomes relevant only after « 10 6 s after the dynamic merger phase. Post-merger evolutions 
have been considered for short and long SOF cases (Tabled and these differ in the amount 
of n-capture elements originating in the SOF (see position of vertical lines in Fig.HJ). Only 
in the long-SOF case do s-process element enhancements originate in the SOF, and provide 
abundances of n-caputure elements that agree with RCB star observations. 

In the high-T SOF case 13 C is destroyed by He-burning within 10 4 seconds, and although 
iV n reaches a higher peak value of ~ 10 14 cm -3 , it is immediately reduced by neutron poisons, 
such as 14 N and 14 (Eq. 11, Fig. 15, paper 1). Later at 10 5 seconds, the neutron density once 
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again increases and reaches a lower peak value of ~ 6.31 x 10 10 cm~ 3 through the neutron- 
source reaction 17 0(a, n) 20 Ne. But these neutrons are also not available for the production 
of s-process elements due the presence of the 14 N neutron poison. As a result s-process 
elements and elements such as Ca and Ti are unaffected in the high-T SOF and remain close 
to their initial values (Fig.HJ). 



3.2. Stellar and nucleosynthesis evolution 

3.2.1. The evolutionary tracks 

The evolutionary tracks of the four cases start at low luminosity in the HR diagram and 
are similar to each other (Fig.[5J). The thin lines in the lower portion of the plot correspond to 
the initial H-shell burning period and are not expected to reproduce all details of the ongoing 
thermal readjustment of the post-merger star. Fortunately, the details of this part of the 
evolution are not important for the RCB star surface abundance predictions. However, the 
thicker lines indicate when He-shell burning takes over, and the star is now in hydrostatic and 
thermal equilibrium and evolves on the nuclear timescale of the He-burning shell. Observed 



( Clavton 


1996; 


Pandev et al. 


2008) 



The models spend between 0.1 (case 3) and 2.75 x 10 5 yr (case 4) in the RCB phase. 
This age is obtained when the mass- loss rate is ~ 10~ 7 M Q /yr. Increasing the mass-loss rate 
decreases the time spent in the RCB phase. A mass-loss rate of ~ 10 -5 M /yr would reduce 
the period of the RC B phase of th e model to ~ 10 4 years, which is the expected period of 



the RCB star phase flClavtonl 120121 ) 



3.2.2. Nucleosynthesis and mixing: origin of 18 and 19 F 

We discuss the different phases of burning and associated nucleosynthesis processes that 
occur in the post-merger evolution of our models by focusing on case 1 (Table[2]). Dots along 
the HRD track in Fig. [5] correspond to the profiles shown in Fig.|6j Figs. [6]^i)a and b show 
the initial abundance profiles of additional isotopes (compared to Fig.[TJ for case 1. The 
second set of panels (corresponding to the cyan dot in Fig. E]) shows CNO-cycle burning of 
the hydrogen contributed to the merged CO+He WD system from the H-rich envelope of the 
original He WD, which was uniformly distributed throughout the He WD in our construction 
of the initial model (cf. Section I2.1.ip . The maximum temperature during this phase is 
~ 10 7 K, and the H-shell leaves 14 N (via 13 C(p, 7 ) 14 N) and 13 C (via 12 C(p, 7 ) 13 N(/3+) 13 C) 
behind. This phase results in an overall increase in the abundances of 13 C and 14 N in the 
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Fig. 5. — The HR diagram of the four considered cases. The thinner portions of the lines 
indicate the evolution when the star burns H in a shell and the grey box shows the area of 
the HR diagram in which RCBs are observed to lie. The dots correspond to timesteps in the 
evolution: (i) - years, empty black circle, (ii) - 4.2 x 10 5 years, cyan dot, (iii) - 6.7 x 10 5 years, 
violet dot, (iv) - 1.3 x 10 6 years, brown dot, and (v) - 7.2 x 10 6 years, black dot. 
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Fig. 6. — The abundance distributions for nuclear species between 1 H and 22 Ne at three 
timesteps of evolution for case 1. The panels correspond to the dots in Fig. [5j The line of 
every nuclear species connects the abundance of that species at every 20 th zone in the model. 
The panels in the left column have log -Dtotar^fkolid cyan line) plotted in them and those in 
the right column have log T (solid orange line). 




Fig. 7. — The abundance distributions for nuclear species between 1 H and 22 Ne at two 
timesteps of evolution for case 1. The panels correspond to the dots in Fig. [5j The line of 
every nuclear species connects the abundance of that species at every 20 th zone in the model. 
The panels in the left column have log Dtotai (solid cyan line) plotted in them and those in 
the right column have log T (solid orange line) . 
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envelope. 

Eventually, the temperature in the star increases further to ~ 10 8 K, and beginning at 
m m i n , He-shell burning starts (violet dot in Fig. [5]). By this time, the H-shell has almost 
completely burnt outwards to retreat to a thin layer close to the surface (Fig. Etli)a). The 
first neutron burst, via a-capture on 13 C is triggered at this time, and consequently the 
13 C-shell moves outwards as it gets burnt to 16 O. 

In the time between panels (ii) and (iii) of Fig. [6] (violet and brown dots in Fig. [5]), partial 
He-burning of 14 N raises the 18 to its peak value, where T ~ 1.20 x 10 8 K, while some of 
this 18 O be gins to get destroyed to 22 Ne by further a-capture closer to 77i m in- 

As He-burning starts, the 13 C(a,n) 16 reaction releases neutrons at a low density of 
N n fa 10 5 cm -3 . This is too little to produce s-process elements, in particular in view of the 
large 14 N abundance that acts as a neutron poison for the s process. However, the neutron 
poison reaction enables the production of 19 F. The abundances of 15 N and 19 F increase 
via the reaction chain 14 N(n, p) 14 C(p, 7) 15 N(a, 7) 19 F. Panel (iii) shows how this reaction 
chain is able to produce 15 N and 19 F efficiently in this simultaneous slow mixing and burning 
environment. The challenge for any scenario for the simultaneous production of 19 F and 18 
is the relative efficiency of the a-capture, which is increasing in the order 13 C, 18 0, 15 N and 
finally 19 F. 15 N appears in the more external layers, where 13 C releases neutrons. However, 
19 F can only form from 15 N in deep enough layers, where typically 18 would already be 
destroyed. Therefore, mixing must be fast enough to bring 15 N below the region where 18 
can form, but slow enough to not transport all that 18 down where it would react with a 
particles into 22 Ne. 

This balance between mixing processes and nucleosynthesis is also realized in the ramp 
region (Eq.|2c]) of our adopted mixing profile. Such a continuous transition from the envelope 
mixing regime to the non-mixed core is a physically reasonable assumption, and it has to be 
emphasized that without such a ramp transition our models could not reproduce the observed 
abundances of F as well as they do. It is this ramp transition region that is responsible for 
much of the continuous increase in 19 F during the later part of the evolution of our model 
stars to the RCB star regime (panel (iv) and (v) in Fig. [7J black dot in Fig. [5]). 

If we assume that the diffusion coefficient drops more steeply at (via the parameter 
/ in Eq.[2c]), then less 19 F is produced. For a larger / parameter, the mixing profile extends 
deeper and the surface amount of 18 decreases, while that of 16 increases. Therefore, 
the combined abundance evolution of 18 and 19 F provide strong limits on the width of the 
transition region between the mixed envelope and unmixed core at and below the He-burning 
shell. With our choice of / the models are allowed to produce the maximum amount of 19 F 
without destroying too much 18 and increasing the 16 at the surface. The ramp-transition 
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feature is thus a sensitive, integral and unavoidable part of the mixing profile. 

In panel (iv), Fig.[7J T has risen at m min to ~ 3.16 x 10 8 K, and the triple-a chain has 
been fully activated causing 12 C to increase and a further a-capture causing 16 to increase 
(Fig-Etlv)a). Below m^, at which the mixing profile drops, complete He-burning of 18 and 
19 F to 22 Ne takes place. By the time the star enters the RCB phase in the HR diagram, it has 
a strong He-shell burning (panel (v) of Fig. black dot in Fig.[S]). Between panels (iv) and 
(v), a second weaker neutron-source, that of 22 Ne, is initiated through 22 Ne(a, n) 25 Mg. The 
neutron densities {N n ) during the post-merger evolution range between 10 4 and 10 6 cm -3 . 

The decrease in mixing efficiency follows the He-burning shell by design (cf. Section [2731 
Eqs.[2]). This design prevents the material that is processed by strong He-shell burning, such 
as 12 C and, particularly, 16 0, from reaching the surface (panels (iv)a and (v)a, Fig.[7|), where 
it would cause surface abundance changes that are not observed. The design also prevents 
18 and 19 F isotopes from above to be mixed to the higher temperatures in the interior, 
where they would be easily destroyed. Thus, the formulation of the mixing profile avoids 
an excess pollution of the surface by 16 and preserves the 18 and 19 F abundances in the 
envelope. As explained in Section I2.3[ the fast decrease in mixing efficiency in our mixing 
model is motivated by the positive entropy gradient induced by the He-burning shell. 

The results of similar simulations for the other three cases from Tabled are shown in 
Figs. and they are compared with observations in Section l3".3.2l 

3.2.3. s process in the post-merger evolution 

As seen in the previous section n-capture nucleosynthesis plays an important role in 
producing the observed F abundance. In case 4 (q = 0.5, high-T SOF, no SOF s process 
prodcution, cf. Section 13. ip the neutron exposure in the post-merger evolution is sufficient 
to produce s-proecess elements. In this lower-g case we have deduced a higher dredge-up 
mass from the dynamic merger simulation which implies a higher initial abundance of 12 C. 
The lower He- WD mass of the low-g case implies a larger WD envelope mass and therefore a 
larger proton abundance in the post-merger envelope. Taken together the low-g post-merger 
evolution features a larger amount of the n-source 13 C. This leads to a higher neutron density 
in case 4 (N n fh 10 6 cm -3 ) compared to cases 1-3 (N n « 10 4 — 10 5 cm -3 ) (Fig. |SJ). Thus, in 
the post-merger evolution of case 4 the abundances of s-process elements are increased by 
up to ~ ldex (e.g. in the case of Y) from their initial values. In cases 1-3 the post-merger 
evolution n-capture nucleosynthesis is significantly less efficient. 
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Fig. 8. — The abundances by number fraction of 13 C/100, Y, Ba and Ti during the post- 
merger evolution of cases 1 and 4, along with neutron densities. The snapshots of time for 
these profiles are at 3.1 x 10 5 yr (case 4) and 4.5 x 10 5 yr (case 1), between the cyan and the 
violet dots in Fig.[5j 



3.3. Comparison of the surface abundances in our models with observations 



RCB stars a re classified into two categories based on their values of [Fe] and Si/Fe 
and S/Fe ratios ( Lambert Sz Rao 1994h . RCB minority stars (plotted as black squares, 
Fig. ED are identified as having lower [Fe] and higher Si/Fe and S/Fe ratios than the majority 
(plotted as black stars , Fig.[9]). We adopt the data of measured elemental abundances of 
RCBs summarized by iJeffery et al.l (120111 ) and use the solar scaled log representation of 
abundances, i.e. [X]= — Cj where q = log rij + C, rii being the relative number abundance 
of species i. This data summary does not include the individual error bars on each data 
point, but in general the errors are quoted to be 0.2-0.3 dex. 



3.3.1. C> N> : constraints on the mixing 

The abundance of C+O traces He-burning, while the C/O ratio signals the relative 
importance of the triple-a to 12 C(a,7) 16 reactions which is sensitive to the temperature. 
The elemental abundance of N is the result of CNO processing. Therefore, the (C+0)/N 
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Fig. 9. — The (C+0)/N ratio of the four computed cases (case 1 (orange symbols), case 2 
(red symbols), case 3 (blue symbols) and case 4 (green symbols)) and in RCBs (minority - 
black squares, majority - black stars) against their C/O ratio. Two minority stars have been 
omitted due to their C/O ratio being out of the range of the plot. R CrB is indicated by 
the magenta star in the plot. Plotted for every case is the surface value of these ratios when 
the model is within the RCB range (empty circles), when it enters the RCB phase in the 
HR diagram (Fig. [5]) (empty triangles) and when it leaves this phase (filled triangles). All 
observed RCBs that have N<0 in their surface are indicated by filled symbols. 
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ratio indicates the relative importance of He- to H-burning processes in determining the 
observed abundance distribution in RCB stars. Their distribution in the (C+0)/N vs. C/0 
plane shows that the majority lie in a particular range of (C+0)/N ratio, of 1-3 and are 
more spread in their C/O ratio, with most RCBs having a value of 1-10 (Fig.[9]). We find 
that on keeping the amount of mixing constant at the initial value of -DLag in our models, the 
14 N in the envelope reduces continuously by being burnt to 18 and this causes an overall 
decrease in the mass fraction of nitrogen at the surface. On continuing the evolution of 
the model at the same rate of mixing up to the end of the RCB phase, the (C+0)/N ratio 
increases at least by up to a factor of 3 higher than its value in most RCB stars (except for 3 
stars in Fig.[9]but these have N>0). If the rate of mixing is lowered and restricted to smaller 
depths from the surface once the model attains the (C+0)/N and C/O ratios of most RCBs 
(the empty circles in Fig.|H]), the 14 N on the surface is preserved through the RCB phase 
of the star. Thus in order to attain the distribution of C, N and O of RCBs, our models 
strongly in dicate that the magn itude of mixing and its reach from the surface must decrease 



with time. ISchwab et al.l (120121 ) find that continuous transportation of angular momentum 
from the merged system causes differential rotation to reduce in the envelope. Although 
the timescale over which they see this happens in their simulations is shorter than what our 
models indicate, the idea that such a reduction in mixing driven by rotation can happen is 
very feasible. 

The difference in the surface abundances between restricting mixing at the time when 
the model lies in the observed domain of (C+0)/N and C/O and on continuing mixing as 
the star enters the RCB phase, is not significant for elements between Na and La. However 
a considerable difference of the order of 0.5 dex does arise in the C, N, O and F abundances 
as seen in in Fig.[10j Case 4 does not at any point of time in its evolution have a model with 
C>N>0 on its surface (hence the absence of circles in Fig. ITU]) and its O-content is almost 
the same as the C-content in its surface. 



3.3.2. Elemental surface abundances 

In this section we provide an element by element discussion based on the following 
diagrams. Fig. [HI] provides a summary for the surface abundance evolution through the 
RCB phase of all four cases. Fig.[TTJ compares the O-isotopic ratio and F abundance of 
models with observations. The [X] values of the observed nuclear species in majority and 
minority RCBs are shown in Fig. [121 along with the lowest and highest values of each element 
from the four cases, between the time when star attains the observed C/N and N/O ratios 
in RCBs and when it exits the RCB phase in Fig.[SJ 

As RCBs are associated with the old bulge population (intermediate Pop II, Sec. CP, we 
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Fig. 10. — Surface [X] abundances of case l,2©ase 2, case 4 and case 4. The empty squares 
stand for the initial envelope abundance while the meaning of the rest of the symbols are as 
explained in Fig. [9j For species between Mg to La, the difference in abundances between the 
models that enter and exit the RCB phase is < 0.2 dex, hence only the average abundances 
of these species in this phase are plotted (filled triangles). 
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Su rface 16 0/ 18 O ratios (squar es are the observed dat a ; emp ty squares from 



mpt 

Clayton et al.l (l2007f ) and filled ones fr om iGarcia-Hernandez et al.l (l2010l )) and [F] (black 



diamonds are the observed values from IPandey et al.l (120081 )) against T e $ for the four cases 
(colours are the same as in Fig. ITUj) . The model abundances of the nuclei plotted here 
(hexagons) are averaged between the timesteps at which the star enters and leaves the RCB 
phase (Fig. |5j). 
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assume that their range of Fe abundances reflects their initial metallicity. Our models were 
constructed with solar metallicity, which allowed us to take advantage of our detailed nucle- 
osynthesis and stellar evolution progenitor models, currently available only at this metallicity. 
Fortunately, this does not affect the relevance of our models for comparison with observa- 
tions in any significant way. The nuclear species in our models can be divided into primary 
and secondary types, depending on the nature of the processes that affect their abundances. 
With this difference in mind our models can be compared with observations of RCB stars 
for a range of metallicities in the following way. 

The primary elements are formed entirely from the burning of H and He and so, are 
independent of the initial [Fe] abundance of the star. Therefore, in Fig.[T2J one must compare 
only between the [X] values of these elements between the models and the observations. 

The abundances of secondary elements in the model depend on the proportion of their 
seed nuclei which were formed in previous generations of stars. These include species between 
Na and Ti and the s-process elements. Therefore, in Fig. [121 one must compare the offset of 
the [X] value of the element with the solar scaled composition shown as a dashed green line. 

H and Li: The high amount of hydrogen depletion reported for all RCBs, is a natural 
outcome of our models. But Li being a fragile element is easily destroyed by H- and He- 
burning during the merging phase of the WDs itself and consequently in our models too. 
Therefore we do not include H and Li in the discussion of surface abundance results. 

Primary elements 

C, N, O and F: One of our main objectives in this work was to find the characteristic high 
enhancements in F and the low 16 0/ 18 ratios of RCBs. Fig.lTTl shows the good agreement 
between our models and the observed range of these species. Three of the cases have 16 0/ 18 
ratios between 6 and 15 that are within the measured range in RCBs, while this value is 
higher by a factor of 2 in case 4. The model F abundances are as well spread over nearly 
the whole range of observed values with the maximum model value of 2.35 dex compared to 
the observed maximum of 2.7dex. The 12 C/ 13 C ratios range from 230 (case 4) to ~ 2500 
(case 3) and are therefore always higher than the observed lower limit of 100. 

The C, N, O and F model surface abundances rely exclusively on the nucleosynthesis and 
mixing during the post-merger evolution phase. Independent of the choice of the progenitor, 
the surface C, N and O values of the four cases are within the same order (Fig. fTOl) . The 
surface F abundances are close to each other between cases 1, 2 and 4 (within a factor of 
1.5, Fig. fTT]) . The highest F abundance is obtained in case 2 since some of it was generated 
in the SOF by He burning, while the lowest is in case 2 due to a much smaller depth of 
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mixing, which led to a smaller mass available to produce it. Case 4 has the highest amount 
of since its initial SOF was burnt at 2.42 x 10 s K, thereby increasing the initial 16 0. This 
is reflected in the high 16 0/ 18 ratio ~ 41 of this well (Fig. llip . In summary, the 

model C, N, O and F abundances lie as well in the observed range. 

Secondary elements 

Na - Ti: Both Na and Mg are enhanced compared to initial surface values in all four cases. 
Na is formed by proton captures on isotopes of Ne during the post-merger evolution. Thus, 
as the H-shell moves outwards in our simulations (Fig. Ep)a), the abundance of Na increases 
throughout the envelope. Na is also initially enriched in case 3 due to its production in the 
A-TP AGB progenitor used to construct the CO WD of this case. 

For the rest of the species from Al to Ti, the abundances are independent of the choice 
of the progenitor AGB model. In cases 1-3, these elements are changed in the low-T SOF, 
depending on its duration of burning, whereas in case 4 their abundances are changed during 
the post-merger evolution phase (as mentioned in section 13.11) . 

Al is first generated by proton-capture on Mg isotopes and later, being a neutron poison, 
is destroyed by (n,p) reactions once again to Mg via, 26 Al(n, p) 26 Mg. The Mg content also 
increases by proton-capture on 23 Na. Thus, although Mg is increased in all 4 cases, the level 
of reduction of Al depends on the interplay between the temperature and burning timescale 
in the SOF, which explains why the Al abundance can be either higher (cases 1 and 2) or 
lower than its initial value in the envelope (cases 3 and 4). 

The abundance of Ca decreases in all cases, due to n-capture reactions on its isotopes. 
40 Ca(n, 7) 41 Ca is a significant n-capture channel and the 41 Ca thus formed is used up by fur- 
ther (n,p) and (n,a) reactions, thus decreasing Ca. The reaction chain 44 Ca(n, 7) 45 Ca(/3~) 45 Sc 
followed by 45 Sc(n, 7) 46 Sc(/3~) 46 Ti further reduces Ca and in the process increases the amount 
of Ti. 

From Fig. [T21 the range of Na values from our models is consistent with the data trend, 
while those of Al are in marginal agreement with the enrichment levels seen in RCBs. Since 
there are only 6 stars which have Mg measured and as they are well scattered, it is difficult 
to draw out a trend in this element; we only barely get the upper limits and do not get the 
decrements, as seen in two of those stars. Upper limits are also closely achieved in Ti and 
Ca, while our lower limits do not reach down to those found in some RCBs. 

Ca is seen to be enriched in a few RCB stars with the lowest [Fe], which none of our 
models show. Likewise, Si and S hardly change in all four temperatures are not 

hot enough to produce them in the initial SOF or in the star itself. Thus our models do 
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not show the overabundances of Si and S observed in RCB stars with the lowest [Fe] (see 
discussion in Section @J. 

s-process elements: The abundances of s-process nuclei depend on the seed Fe nuclei 
available in the star. Although the neutron source is primary, the main neutron poison 
is also primary, and therefore the s process is secondary. Thus, one would expect [X] of 
s-process elements to rise rise more or less in sync with [Fe]. 

Among our simulations the low-g case with a high-T SOF (case 4) shows enhanced 
abundances of s-process elements from the neutron flux available during the post-merger 
evolution phase (Section 13. ip . while the high-g case with low-T and long SOF (case 2) 
produce the s-process elements in the SOF before the post-merger evolution (Section |3.2.3[) . 
Cases 1 and 3 produce only small amounts of s-process elements. 

With few exceptions almost all RCB stars show some enhancement in s-process elements 
on their surface, compared to the solar-scaled abundance level. The observed range of 
enhancements agrees with the range of our four cases, except for a few stars that seem to be 
depleted in s-process elements. 

4. Discussions and Conclusions 

Our post-merger models are based on realistic progenitor evolutions with complete nu- 
cleosynthesis predictions and the dynamic merger history presented in paper I. These models 
adopt an envelope mixing model that is motivated by the rotation of the merged object. With 
these assumptions the models reproduce almost all the observed characteristics of RCB stars, 
in particular the anomalous isotopic signatures of O, as well as C and the F enhancements. 

Model predictions are also in good agreement with the observed abundances of s-process 
elements. These do not originate in the progenitor AGB star. Depending on the mass ratio 
of the merging WDs the s process originates either in the SOF, which is a characteristic of 
the dynamic phase (case 2) or in the post-merger phase (case 4). In the former case s-process 
enhancements are only seen if the SOF is assumed to persist for more than 10 6 s (long SOF). 

The enhancements of C, N, O and F are due to nucleosynthesis in the post-merger 
evolution of the star. The initial abundances constructed from the progenitor AGB and 
RGB models are important with respect to the amount of H available, and the profile of C 
and O into which dynamic and post-merger dredge-up mixing may reach. Despite adopting 
a range of progenitor AGB models, we find that the C, N, and O abundances are similar for 
all cases. 
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Our models also lie within the observed range of luminosity and effective temperature 
of RCB stars. The differences between the surface abundances when the models enter and 
leave the RCB phase are < 0.5dex for C, N, O, and F, and < 0.2dex for the other nuclear 
species (this is less than the observational errors). Thus, our models can also account for 
the four hot RCB stars between 15000 and 20000K. The models show a range of 16 0/ 18 
and F in agreement with observations. The 12 C/ 13 C ratios of the models are higher than the 
observed lower limit for RCB stars. 

The model parameters are constrained by the diverse range of observational properties. 
In addition to the depth of mixing during He-shell burning phase of the post-merger evolution 
(k), the mixing limit m min depends on the AGB progenitor. For an advanced AGB progenitor 
model the depth of mixing must be limited to avoid an excessive 16 dredge-up to the surface 
since the outer layers of the evolving AGB CO-core get increasingly O-rich. Other parameters 
of the mixing model (Section 12. 3j) calibrated to reproduce the observed surface abundances 
are the width of the partially mixed region below the proper mixing regime in the envelope 
and the efficiency of mixing. The former needs to be large enough to produce enough F but 
small enough to avoid destruction of 18 O. The latter must be large enough to mix enough 
18 to the surface but small enough to avoid contamination with 16 0. Surface abundances 
are also determined by the SOF properties, which are mostly a result of the mass ratio of the 
merging WDs. The duration of conditions for nucleosynthesis in the SOF is uncertain and 
has been approximated by assuming two long and short SOF. Finally, our models 

suggest that the adopted additional envelope mixing should subside just before the RCB 
star phase is reached. Without such a decrease of mixing efficiency with time, continuous 
mixing throughout the evolution of the star causes the star to become O-rich as against 
being C-rich. 

Although in the end the choice of parameters (progenitor model, SOF assumptions, 
mixing model) appears to have a lot of parameters, it turns out that all of these seem to 
be identifiable with uniques nucleosynthetic and mixing signatures. Due to large number 
of nucleosynthesis processes, each with their individual dependence on temperature and 
on other, related nucleosynthesis processes (as for example the neutron production and 
consumption for the F production), most of these parameters are therefore well constrained 
in our models. 

A legitimate question would be how degenerate our solutions are. In view of the many 
cases that we have investigated, many of them not described in detail in this paper, we are 
not aware of any alternative way to arrive at such a level of agreement between observations 
and models. However, one has to admit that significant uncertainties, especially regarding 
the properties and evolution of the SOF as well as non-spherically symmetric aspects of the 
post-merger mixing model remain. 
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Observations of RCB stars show some diversity and our four cases reproduce a spread 
in abundance features. Most of this spread in model predictions is in fact due to differences 
in the cases, such as the mass ratio of the WDs, or the progenitor model, which represent a 
real source of diversity corresponding to the range of WD mergers that are expected to end 
up as RCB or HdC stars. 

When comparing obser yations and model prediction it must be kept in mind that dust 
depletion may play a role (lAsplund I Il998l : lAsplund et al.l l2000l ). After all the irregular 
variability of RCB stars is associated with dust enshrouding. This possibility cannot be 
ruled out, in particular to explain the distinctly low [Fe] values of the minority RCBs from 
the majority. It does not appear that the dust-gas separation in RCB stars. An ISM-like 
dust depletion can explain the high [S/Fe] and [Na/Fe] observed in RCBs, but not their 
[Al/Fe] and [Ca/Fe] values, both of which are much lower in the ISM. In particular, the 
Si in a carbon-rich dust environment is expected to condense into SiC grains, thus severely 
reducing the Si content. However all RCBs have high [Si/Fe] ratios. It should be noted 
however, that only amorphous graphite grains and no SiC grains have been detected in 
RCBs. This maybe indicative of a different dust-gas separation mechanism from that of the 
ISM. 

An alternative interpretation of the large Si and S overabundances seen in the RCB 
stars with the lowest [Fe] considers the mass of the merging WDs. Si and S are produced in 
O bur ning at T ^ 10 9 K, which could be achieved in the SOF of a merger of more massive 



WDs (lLoren-Aguilar et al.ll2009L T — 0.6 — 1.0 x 10 9 K for combined mass > 1.2M ). The 
theoretical initial-final mass relationship predicts, for a given initial mass, higher WD masses 
for lower metallicity. Assuming the initial mass function does not depend on metallicity in 
the range relevant for RCB star progenitors, higher WD masses are expected for lower 
metallicity mergers, and this may be the reason for the larger enhancements in S and Si for 
the RCB stars with the lowest Fe content. Such higher- mass RCB stars would naturally 
also be more luminous, and this hypothesis therefore predicts that RCB star luminosity is 
correlated with Si and S enhancements. 

On a preliminary note, we do see enhancements of 22 Ne on the surface of our models 
that seem to agree with the levels of Ne found in EHes. A detailed comparison of our models 
with the observations of EHe is, however, beyond the scope of this work. 

Finally, a few exceptions to the usual RCB star properties inclu de those with higher 



than average H, in cluding four RCB stars in which Li is enhanced (lAsplund et al.l 12000 



Kipper et al.l 120061 ) , as well as four other RCBs that have a re latively higher amount of 
13 C at their surface faao fe Lambert II2OO8I : lAsplund et allboooh . with 12 C/ 13 C< 25. Nei- 
ther of these p r opert ies are not found in our models. Recent SPH simulation results of 
Longland et al.l (120121 ) indicate that, to be preserved, any Li that maybe produced in a DD 



36 



merger will have to be ejected into the mass that forms clouds around the RCB star. Since 
mass-ejection during the dynamic merger phase is not part of the hydrodynamic simulations 
of paper I, our models cannot predict lithium abundances from such a scenario. Thus, the 
question of whether Li-rich RCBs are an outcome of the FF scenario or of the merger between 
WDs remains open and needs a separate, more in-depth analysis. 
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